Electromagnetic Decay of the S°(1385) to A7 
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The electromagnetic decay E°(1385) —5- A7 was studied using the CLAS detector at the Thomas 
Jefferson National Accelerator Facility. A real photon beam with a maximum energy of 3.8 GeV 
was incident on a proton target, producing an exclusive final state of K + T,*°. We report the 
decay widths ratio E°(1385) -> A7/ E°(1385) -)■ Att° = 1.42 ± 0.12(stat)t^7(sys)%. This ratio 
is larger than most theoretical predictions by factors ranging from 1.5-3, but is consistent with the 
only other experimental measurement. From the reported ratio we calculate the partial width and 
electromagnetic transition magnetic moment for E°(1385) — > A7. 



PACS numbers: 13.40. Em,14.20. Jn,13.30.Ce,13.40.Hq 
INTRODUCTION 

One well-known success of the constituent quark model 
(CQM) is its prediction of the low-mass baryon magnetic 
moments, using just the SU(6) wavefunctions [l|,|2|. Cal- 
culations of the magnetic moments assuming that 
quarks behave as point-like Dirac dipolcs, are typically 
within ~10% of the current measured values [4j. How- 
ever, today we know that the spin of the proton is much 
more complex than the CQM representation, with only 
about one-third of the proton's spin coming from the 
quarks and the rest of the spin resulting from a combina- 
tion of the gluon spins and the orbital angular momen- 
tum of the quarks 5, 6]. Clearly, the CQM is an over- 
simplification of the spin dynamics inside baryons, yet 
somehow the CQM captures the degrees of freedom that 
are relevant to the measured magnetic moments. Further 
measurements of baryon magnetic moments, utilizing the 
electromagnetic decays of excited baryons, will continue 
to test our understanding of baryon wavefunctions. 

Experimentally, it is difficult to measure the electro- 
magnetic (EM) transitions of decuplet-to-octet baryons 
because of the competition between EM and strong de- 
cays. For example, the branching ratio for EM decay of 
the A resonance has been measured to be about 0.55% 
[ij and branching ratios for other decuplet baryons are 
predicted to be of the same order of magnitude. The 
EM transition form factors for the A may be directly 
measured via pion photoproduction 0, [1[ . 

It has been shown Q that pion cloud effects contribute 
significantly (~40%) to the jp — » A + magnetic dipole 
transition form factor, Gm(Q 2 ), at low Q 2 (below ~ 0.1 
GeV 2 ). In the naive non-relativistic quark model [ldj . 
the value of Gm(0) is directly proportional to the pro- 
ton magnetic moment, and measurements of Gm near 
Q 2 = can only be explained (within this quark model) 
if the experimental magnetic moment is lowered by about 
30%. This again suggests that the CQM is an over- 
simplification of reality. 



To extend these measurements to the other decuplet 
baryons, which have non-zero strangeness, hyperons must 
be produced through strangeness-conserving reactions. 
Then their EM decay, which has a small branching ra- 
tio, must be measured directly. Although these measure- 
ments are difficult, it is important to measure the EM 
decays of strange baryons to extract information on their 
wavefunctions, which in turn, constrains theoretical mod- 
els of baryon structure. The measurements of EM tran- 
sition form factors for decuplet baryons with strangeness 
may also be sensitive to meson cloud effects. A com- 
parison of the EM decay measurements to predictions of 
quark models for decay of decuplet hyperons, S*, to octet 
hyperons, Y , can provide a measure of the importance of 
meson cloud diagrams in the E* — > Yj transition. 

Here, we present measurements of the EM decay 
E*° — > A7 normalized to the strong decay S*° — > Air . 
The present results can be compared to previous mea- 
surements of the E*° EM decay [ll| (also from CLAS 
data) that had a larger uncertainty (^25% statistical 
and ~15% systematic uncertainty). The smaller uncer- 
tainties here are due to a larger data set (more than 10 
times bigger) and, subsequently, better control over sys- 
tematic uncertainties. The reduced uncertainty is impor- 
tant because, as mentioned above, meson cloud effects 
are predicted to be on the order of ~30-40%. In order 
to quantify the effect of meson clouds for baryons with 
non-zero strangeness, it is desirable to keep measurement 
uncertainties below ~K)%. 

There are many theoretical calculations of the EM de- 
cays of decuplet hyperons such as: the non-relativistic 
quark model (NRQM) IT2I [l3| . a relativized constituent 
quark model (RCQM) [14J, a chiral constituent quark 
model (xCQM) 0, the MIT bag model 0, the bound- 
state soliton model [l7j], a three-flavor generalization 
of the Skyrme model that uses the collective approach 
, an algebraic model of hadron structure [20j , and 
heavy baryon chiral perturbation theory (HB%PT) [2ll |. 
among others. Table |T] summarizes the theoretical pre- 
dictions and experimental branching ratios for the EM 
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TABLE I: Theoretical predictions for the models referenced in the text and the experimental values for the electromagnetic 
decay widths (in keV). 



Model 



A(1232) -> N-y 



£(1193) -» A7 



NRQM [12, 13,^6] 
RCQM [14] ~ 
X CQM [15] 
MIT Bag [16J 
Soliton [17] 
Skyrme [18, ,19] 
Algebraic model [20] 
HBxPT [21] t 



360 
350 



309-326 

341.5 
(670-790) 



Experiment [4J 



8.6 
4.1 

4.6 



8.6 



£(1385) ->■ A7 



273 

267 

265 

152 

243 
157-209 
221.3 
290-470 



660±47 



T Normalized to experiment for the A — > N'y range shown. 



8.9±0.9 



470±160 



transitions of interest. 



THE EXPERIMENT 



A comprehensive study of electromagnetic strangeness 
production has been undertaken using the CLAS detec- 
tor at the Thomas Jefferson National Accelerator Facil- 
ity. Many data on ground-state hypero n photoproduc- 
tion have already been published 2214241] using data from 
the CLAS run group gl and all data sets. The gl exper- 
iment had an open trigger |22j and a lower data acqui- 
sition speed, whereas the gll experiment required that 
at least two particles be detected 24 1, and used a higher 
beam current, resulting in much higher data acquisition 
speed. The result is that the gll data set had over 20 
times more reconstructed events than the gl data set. 
The present results use the gll data set, whereas Taylor 
et al. I 1 11 used the gl data set. Published CLAS re- 
sults |24j from gll demonstrate the accurate calibration 
of this data set, and that the cross section of jp — > K + A 
matches previous CLAS data. 



The EM decay of the E*° is only about 1% of the total 
decay width. To isolate this signal from the dominant 
strong decay S*° — > A71- , the missing mass of the de- 
tected particles, jp — > K + A(X), is calculated. Because 
of its proximity to the it peak in the mass spectrum from 
strong decay, the EM decay signal is difficult to separate 
using simple peak-fitting methods. The strategy here is 
to understand and eliminate as much background as pos- 
sible using standard kinematic cuts, and then use a kine- 
matic fitting procedure for each channel. As described 
below, by varying the cut points on the confidence lev- 
els of each kinematic fit, the systematic uncertainty as- 
sociated with the extracted ratio for EM decay can be 
quantitatively determined. The increased statistics for 
the gll data helps greatly to study and to determine the 
systematic uncertainty associated with the measurement. 



For the present measurements, a bremsstrahlung pho- 
ton beam was produced from a 4.019 GeV electron beam, 
resulting in a photon energy range of 1.6-3.8 GeV. The 
photon energy was deduced from a magnetic spectrom- 
eter [25| that "tagged" the electron with an energy res- 
olution of ~ 0.1%. A liquid- hydrogen target was used 
that was 40 cm long and placed such that the center of 
the target was 10 cm upstream from the center of CLAS. 
As mentioned above, a trigger requiring two charged par- 
ticles in coincidence with the tagged electron was used. 
The data acquisition recorded approximately 20 billion 
events. Details of the experimental setup are given else- 
where Hm. 



Event Selection 

We selected events for the reaction -fp — > K + Y,*° , 
where the S*° decays with 87. Oil. 5% probability to A7r° 
and 1.3±0.4% probability to A7 [4]. The A then decays 
weakly with 63.9±0.5% probability to pn~ leading 
to the final states jp — > K + pir~ir° and 7p — > K^pn^j, 
respectively. The charged particles are tracked by the 
CLAS drift chambers through the magnetic field of the 
spectrometer, giving their momentum, and are detected 
by the time-of-flight scintillators, giving their velocity. 
The drift chamber tracking covariance matrix is obtained 
for each track. This contains the uncertainty in each mea- 
sured variable used in track reconstruction along with the 
appropriate correlations. The ir° and 7 must be deduced 
indirectly using conservation of energy and momentum 
via the missing mass technique. 

In the present analysis, two positively charged parti- 
cles and one negatively charge particle are selected. The 
mass of the detected particles was calculated from the 
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1: Missing mass squared (Mjj) for the reaction -yp — > 
p-K + -K~ (X), where the ir + was a potentially misidentificd 
kaon. Events above the dotted line at 0.01 GeV 2 were kept. 
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measured velocity and momentum. The mass is given by 



FIG. 2: The distance of closest approach between the proton 
and n~ shown in cm. The cut used at 5 cm is indicated by 
the dotted line. 



TOcal 



/p 2 (l-/3 2 ) 



ft 2 



(1) 



where /3 = L/t mcas for path length L and measured timc- 
of- flight tineas, and the speed of light is set to 1. The 
pions, kaons, and protons were identified using mass cuts 
of 0.0 < M„- < 0.3 GeV, 0.3 < M K + < 0.8 GeV, and 
0.8 < M p < 1.2 GeV, respectively. From this initial 
identification it is possible to incorporate additional tim- 
ing information to improve event selection. The time- 
of-flight t moas is the time difference between the event 
vertex time and the time at which the particle strikes 
the time-of-fiight scintillator wall at the outside of the 
CLAS detector. We define At = t mcas — t ca i, where t ca \ 
is the time-of-fiight calculated for an assumed mass such 
that 



^cal 



L\ 1 + 



(2) 



where m is the assumed mass for the particle of interest 
and p is the momentum magnitude. A cut on At or m ca \ 
should be effectively equivalent. 

Using At for each particle it is possible to reject events 
that are not associated with the correct RF beam bunch, 
which are separated by 2 ns. This is done by accepting 
only events with |At| <1 ns. 

A A/3 cut also helps to clean up the identification 
scheme. A/3 is the difference between the above measured 
/3 and the calculated /3 C defined by (3 C = pj \Jp 2 + m 2 , 
where p is the particle momentum and m is the known 
particle mass. The good events were required to have 
-0.02 < A/3 < 0.02. 

The energy lost by charged particles passing through 
the CLAS detector was accounted for by adjusting the 
measured particle's energy according to the average 
dE/dx losses in the target material, target wall, target 
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FIG. 3: The invariant mass of the p-n~ showing the A peak 
with a Gaussian fit giving a a — 1.3 MeV. A cut of ±0.005 
GeV around the PDG mass of the A is used indicated by the 
dotted lines. 



scattering chamber, and the start counter scintillators 
surrounding the target. After correcting for energy loss, 
several kinematic cuts are applied as described below. 

Due to the finite resolution of the measured velocity 
and momentum, in addition to particle decay-in-flight, 
it is possible that some pions could be misidcntified as 
kaons. To clean up the kaon signal for the analysis, 
it is common to recalculate the energy of the identified 
kaon using the mass of the pion. Then the missing mass 
squared is studied for the reaction jp — > pir + ir~ (A), 
where the ir + is actually identified by the above mass 
cuts as a K + . A spike at zero mass squared indicates 
that the reaction jp — > pir + Tr~ is prominent. Most 
particle misidentified events can be removed by cutting 
slightly above zero, as shown in Fig. [I] The events 
above 0.01 GeV 2 are kept as a cleaner sample of the 
7P — > pK + n~ (X) events. Reactions involving decays 
such as p — > ir + ir~ , where the ir + is mistakenly identified 
as a K + , are vastly reduced by this cut. 
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events passing the cut on the A mass. The dotted lines show 
the 1.34 GeV to 1.43 GeV cut used to select the £(1385). 

FIG. 6: Missing mass squared for the reaction 7p 
K + pir~ (X) after all kinematic cuts. 




M x (pjO (GeV) 

FIG. 5: Missing mass for the reaction 7p — > A(X) for events 
passing cuts on the A and £* masses. The dotted line at 
0.55 GeV shows the cut used to remove the 7p — > K + A chan- 
nel. (A looser timing cut was used to illustrate that these 
accidentals are at the K + mass.) 



The four-momentum of the detected A was recon- 
structed from the proton and n~ four-momenta. The 
distance of closest approach (DOCA) from the proton 
and 7r~ four-momenta is found and restricted to be less 
than 5 cm (see Fig. [2]). A Gaussian fit to the pn~ invari- 
ant mass peak shown in Fig. |3]resulted in a a — 1.3 MeV, 
which is consistent with the instrumental resolution. Af- 
ter restricting the A mass to be 1.1157 ± 0.005 GeV, the 
remaining events were used to construct the missing mass 
off the K + , giving the excited-state hyperon mass spec- 
trum shown in Fig. 01 

After making a cut on the E* peak from 1.34-1.43 GeV, 
as shown in Fig. 2J one can study the missing mass off 
of the A, such that jp A(X), shown in Fig. O Small 
peaks are seen at the mass of the kaon and the K*(892). 
The kaon peak is from exclusive 7p — > K + A production 
due to accidentals under the TOF peak, and can easily 
be cut out. The dotted line shows the M X (A) > 0.55 



GeV event selection used to eliminate this background. 

After including all of the cuts listed above, the missing 
mass of the reaction jp — >• K + A(X) is shown in Fig. [5] 
A very prominent peak is seen at the mass of the ir° with 
a very small number of counts about zero missing mass 
due to the EM decay. The counts above the ir° peak 
are mostly due to the jp — > K + YP(X) reaction from 
photoproduction of higher-mass hyperons. 

A small fraction of the events near zero missing mass in 
the spectrum of Fig. [5] come from accidentals and double 
bremsstrahlung. In the case of double bremsstrahlung it 
is possible for false EM decay signals caused by the reac- 
tion 7i + 72P — > K + A + 71 to mimic the final state of in- 
terest 7p — ¥ K + A^f. The 71 from double bremsstrahlung 
will point down the z-axis (along the beam), which can 
also occur if the event is accidental or due to ineffi- 
ciencies in the tagger plane from incorrect electron se- 
lection. By calculating the transverse missing momen- 
tum (P^y — P% + Py), it is possible to eliminate double 
bremsstrahlung. The peak at small values in the distri- 
bution in Fig. [7] was removed by requiring P xy > 0.0009 
GeV 2 as illustrated by the dashed line. Clearly the effect 
is quite small, however this step is critical for an accurate 
measure. 

In Fig. [S] the tail of the 7r° peak continues into the 
zero missing mass region. Resolving the two contribu- 
tions with a simple Gaussian fit does not give systemati- 
cally consistent results. A technique involving kinematic 
fitting is required to separate the background from the 
events of interest, as well as to separate the ir° events 
from the radiative signal. 



SIMULATIONS 

A Monte Carlo simulation of the CLAS detector was 
performed using GEANT [28f . set up for the gll run 
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FIG. 7: (A) Transverse missing momentum and (B) trans- 
verse missing momentum squared for the reaction 7p — > 
K + A(X). The dashed line shows the cut used at Pl y > 0.0009 
GeV 2 . 

conditions. Events were generated for the radiative 
channel (S°(1385) — >• A7), the normalization reaction 
(E°(1385) — > A7T°), and several background reactions; see 
Table HI] for a complete list. Using the data as a guide, 
the photon beam energy dependence of K + production 
and the K + angular dependence were used iteratively to 
tune the Monte Carlo to match the data. After recon- 
struction, the Monte Carlo momentum distributions for 
the proton, ir~ , and K + matched (within error bars) to 
that of the data. The generated Monte Carlo events were 
analyzed using the same analysis procedure used for the 
data. 

After studying the various channels of interest and 
background, a constant t-slope (da/dVl ~ e _bt ) of b =2.0 
GeV -2 was used for the generated — s- if + A(1405) 
channel. The form of the angular distributions of the 
cross section from data were used in the generator to 
produce all the S* simulations. 

KINEMATIC FITTING 

The kinematic fitting employed in this analysis tech- 
nique takes advantage of the information in the measured 
kinematic variables and their uncertainties to fit con- 
straints of energy and momentum conservation, thereby 
improving the measured quantities using constraint equa- 
tions. This procedure is useful to improve the separation 
of signal from background. The method of Lagrange mul- 
tipliers is the approach implemented here to fit the con- 
straints with a least squares criteria 27}. 

Assume there are n independently measured data val- 
ues y, which in turn are functions of m unknown vari- 
ables <7i, with m < n. The condition that y = 
is introduced, where fk is a function dependent on the 
data points that are being tested for each k independent 
variable at each point. 

Because each y^ is a measurement with a correspond- 
ing standard deviation <jfc, the equation y^ = fk{qi) can- 



not be satisfied exactly for m < n. It is possible to 
require that the relationship be satisfied by defining the 
X 2 relation such that 

x2 = E (y fc -/*(< (3) 
fe a k 

and require that the preserved values are qi, which are 
the values of q that minimize x 2 ■ 

The unknowns are divided into a set of measured vari- 
ables (77), such as the measured momentum components, 
and unmeasured variables (u), such as the missing mo- 
mentum or the four-vector of an undetected particle in 
the reaction. The variable Ci is introduced to be used 
for each constraint equation. The Lagrange multipliers 
(Ci) are used to write the equation for x 2 f° r a se t of 
constraint equations T such that, 

X 2 ^ u, C) = (rf - rffV-'irJo - ff) + 2£ T T(rf 1 u), (4) 

where 770 is a vector of initial measured quantities and 
V^ 1 is the inverse of the covariance matrix containing all 
of the resolutions and correlations of the measured vari- 
ables from the drift chamber tracking for each charged 
particle. 

The x 2 minimization occurs by differentiating x 2 with 
respect to each of the variables, while linearizing the con- 
straint equations and obtaining improved measured val- 
ues from the fit. These values are used as the input for a 
series of iterations. The iteration procedure is continued 
until the difference in magnitude between the current x 2 
and the previous value is smaller than Ax 2 est (<0.001). 

The implemented covariance matrix V was corrected 
for multiple scattering and energy loss in the target cell, 
the scattering chamber, and the start counter. These 
corrections to the diagonal terms in the covariance ma- 
trix are applied according to the distance each charged 
particle travels through the corresponding material. 

ANALYSIS PROCEDURE 

A useful kinematic fit for a topology that has a par- 
ticle that is not detected is the 1C fit. Such a fit re- 
quires that a missing mass hypothesis be used to con- 
strain the detected four-momentum, leading to three un- 
knowns from the non-detected particle momentum and 
four constraints from conservation of energy and momen- 
tum. To ensure only high quality A events, an additional 
constraint can be implemented on the proton and ir~ 
tracks to constrain the invariant mass to be the known 
mass of the A. After the detected particle tracks are kine- 
matically fit, the events can be filtered with a confidence 
level cut. In this fit there are three unknowns (p x ) and 
five constraint equations, four from conservation of en- 
ergy and momentum and the additional invariant mass 



7 



condition. No additional constraints are required. This 
makes it a 2C kinematic fit. 

To separate the contributions of the E*° EM decay and 
the strong decay, the events were fit using the hypotheses 
for each topology with the constraint equations, 



T = 



M\ 



-Bboam + M p - E K -E p -E n -Ex 

Pbeam - PK ~ Pp ~ Pit ~ PX 



0, (5) 



where px and Ex are the momentum and energy of the 
undetected 7r° or 7. 

To test the functionality of the kinematic fit used to 
separate the radiative signal from the overwhelming 7r° 
background, the probability density function is used 
to fit the resulting x 2 distribution. The additional con- 
straint on the invariant mass of the A takes the proba- 
bility density function from the more difficult to fit one 
degree of freedom x 2 distribution (containing a singular- 
ity) to the more manageable two degrees of freedom. The 
fit function takes the form, 



fix 2 



Pi 



(6) 



where P2 is a background term, Pi is a quantitative close- 
ness parameter (which gives a measure of how close the 
distribution in the histogram is to the ideal theoretical 
X 2 distribution) , and p) is for normalization. For a kine- 
matic fit to a missing 7 with significant background con- 
tamination from the 7r°, the x 2 distribution will be highly 
distorted. The ideal Pi from a fit to a x 2 distribution 
with no background is determined from simulations. The 
deviation of the Pi fit parameter from the ideal Pi is used 
as an indicator of the signal to background contribution 
going into the kinematic fit under the radiative hypoth- 
esis and how effective a confidence level cut is expected 
to be for that given deviation. 

Using the 7r°-hypothesis for the kinematic fit, the x 2 
distribution follows the trend of the probability density 
function for two degrees of freedom from Eq. ([6]), see Fig. 
[8]A.. The confidence level in Fig. [8j3 is reasonably flat for 
the vast majority of events. The spike at zero confidence 
level in Fig. [8)3 is from events that do not satisfy the 
hypothesis in the kinematic fit. 

For the 7-hypothesis, without any cuts to reduce the 
7T° background, the % 2 distribution is not consistent with 
the expected probability density function for a 2C fit. 
Simulations indicate that the Pi parameter should be 
Pi ~ 0.9, but due to the sensitive nature of the % 2 dis- 
tribution for two degrees of freedom, a fit to obtain the 
Pi parameter does not return a realistic value. This can 
be seen in the distorted shape of the distribution in Fig. 
[9K. Additionally, the confidence level distribution rises 
up near the low confidence end (Fig. [9)3) and is clearly 
not as flat as the distribution in Fig. [8)3. This is an indi- 
cation that the vast majority of data being kinematically 




2 4 6 



FIG. 8: (A) x 2 distribution and (B) confidence level distri- 
bution for a missing it hypothesis in the kinematic fit. 



fit at this stage are not satisfying the base assumption of 
a massless missing particle. This suggests that, even with 
a high confidence level cut, there is still an overwhelming 
amount of 7r° events leaking through. However, it is pos- 
sible to take an additional step in the kinematic fitting 
procedure for a cleaner separation. 

A two-step kinematic fitting procedure is used to sys- 
tematically reduce the large ir° background and optimize 
the extraction of the number of radiative events. First, 
a fit to a 7r°-hypothesis is done and only the low confi- 
dence level {P%(x 2 )) events are retained, followed by a fit 
of these candidate events to a 7-hypothesis and retaining 
the high confidence level (P^(x 2 )) events. Because of the 
previous kinematic cuts, there should now be primarily a 
7T° background and the true EM decay signal. Any other 
background is expected to be very small relative to the 
radiative signal and will be accounted for through simu- 
lations. By first fitting to a 7r°-hypothesis and taking the 
low confidence level candidates, one reduces the proba- 
bility that the surviving candidates will have a missing 
mass of the ir° before they are fit to a 7- hypothesis. 

The selection of the confidence level cuts P^(x 2 ) and 
P^(x 2 ) is derived using simulations. After testing the 
ability to recover various mixed ratios on the order of 
the expected experimental ratio (~ 1%), a Monte Carlo 
(MC) simulation of the data was studied for a given ra- 
tio of the 7p -)• K+Y,* -> K+Air° and jp -> K+Z*° ->• 
A' + A7 channels. The optimization occurs when consid- 
ering both the increase in statistical uncertainty from a 
higher P^?(x 2 ) cut and the increase in MC ratio "recov- 
ery" uncertainty from a lower P£(x 2 ) cut. The "recov- 
ery" uncertainty is defined by the difference in the MC 
generated ratio and the measured ratio or "recovered" 
ratio found by analyzing the MC with a given P"(x 2 ) 
and P^x 2 )- The final confidence level cut in P?(x 2 ) is 
determined by the fit parameter Pi indicating how much 
7T° background is left after the P%(x 2 ) cut. Again, the 
statistical uncertainty and the MC ratio "recovery" un- 
certainty are considered in the optimization of the P^(% 2 ) 
cut. 
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P(x 2 ) 



FIG. 9: (A) The x 2 distribution and (B) the confidence level 
distribution for a missing 7 hypothesis in the kinematic fit 
before the two-step kinematic fit. (C) The \ 2 distribution 
and the (D) the confidence level distribution for a missing 7 
hypothesis in the kinematic fit after the P%(x 2 ) < 1% cu t- 



The results of the optimization study indicate that a 
confidence level cut of P%(x 2 ) < 1% sufficiently reduces 
the 7T° background so that a P 7 (x 2 ) > 10% cut can be 
used to isolate the radiative signal in the kinematic fit to 
7- 

After the two-step kinematic fitting procedure, one can 
again study the 7- hypothesis \ 2 fit- It now looks more 
like a standard distribution for two degrees of freedom, 
returning a value of P x = 0.87 ± 0.06, see Fig. The 
confidence level now appears relatively flat in Fig. [9)3, 
as it should. This is an indication that an improvement 
has been made on the quality of candidates going into 
the fit with respect to the hypothesis. This gives some 
assurance that the candidates going into the secondary 
fit can be accurately filtered with a confidence level cut. 

To ensure the quality of the 7r° extraction, the same 
two-step kinematic fitting procedure is done by first fit- 
ting to a 7 hypothesis and taking the low confidence level 
Pj (X 2 ) candidates, then fitting to the n° hypothesis and 
taking only the high confidence level P%(x 2 ) candidates. 

Once the confidence level cuts are optimized for ex- 
tracting both the 7T° and radiative signal, the final se- 
lected candidates for each case can be seen in the miss- 
ing mass spectrum, see Fig. [TUJ The extracted counts 
are shown for (A) the 7r°, (B) the electromagnetic sig- 
nal, and (C) together in the full spectrum of the missing 
mass squared. The final raw yields taken directly from 
the kinematic fit are n 7 = 635 and n T = 13950. 

The 7T° leakage into the 7 channel is the dominant cor- 
rection to the branching ratio. The final result also needs 
to be corrected for backgrounds, such as K* — > K + X 
and decays to £ + 7r _ , as well as the contributions from 
A(1405) — > Aj. Taking these backgrounds into consid- 
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FIG. 10: 



(A) The n n counts extracted using the confidence 



level cuts P^ < 1% and P° > 10%. (B) The n 7 counts 
extracted using the confidence level cuts P% < 1% and P 7 > 
10%. (C) The counts n v and n 7 shown in the spectrum before 
any kinematic fit. 



eration, and following the notation of Taylor et al. [ll| . 
the branching ratio R = N(Aj)/N(Aw) is 



R 



1 



An^(A 7 )-An 7 AS(A 7 ) 
An 7 ^(AtO + ^L^Ett; 

-An, (^(Att) + ^^(Ett) 



(7) 



where terms starting with A are acceptance factors (given 
below) and 

An, = n OT - N n (A* ->• E+tt") - N W (A* -> £ tt ) 
-N„(A* -> £° 7 ) - AT, (A* -> A 7 ) 
-N n (K* -> Kit ), (8) 



An 7 = n 7 - N 7 (A* -> E+tt") -iV 7 (A* 



E°7r°) 



~A^ 7 (A* -»• E° 7 ) - A^ 7 (A* -> A 7 ) 

-N 7 (K* if 7), (9) 

with n 7 (n„) equal to the yield of the kinematic fits, 
representing the measured number of photon (pion) can- 
didates. In the notation used, lower case n represents 
an observed number of counts, while upper case N rep- 
resents the acceptance corrected or derived quantities. 
Here the it and 7 subscripts indicate the kinematic fit 
hypothesis and the decay channel is shown in parenthe- 
ses (note that A* denotes the A(1405)). These corrections 
are necessary to take into account due to the fact that the 
background underneath the £(1385) is not zero, which 
could lead to an over-counting of the £(1385) contribu- 
tion. For the detector acceptance, the notation has the 
pion (photon) hypothesis from the £(1385) decay given 
by A^ (A^), so that ^(Att) denotes the relative leakage 
of the £ — > Ait decay channel into the A7 extraction 
and A^ (A7) denotes the relative leakage of the A7 decay 
channel into the A-7T extraction. The form of the ratio 
given in Eq. (7) is developed in more detail in the Ap- 
pendix. 



9 



TABLE II: Acceptances (in units of 10 3 ) for the channels used in the calculation of the branching ratio. Here there is a 
P a (x 2 ) < 1% confidence level used with a P h (x 2 ) > 10% cuts. The uncertainties listed are statistical only. 



Reaction 










A(1405) - 


-> 2j 7T 


0.0495±0.0031 


0.001±0.0001 


1.189±0.019 


A(1405) - 


-> £ + 7T- 


0.029±0.002 


0.0013±0.0001 


0.0078±0.001 


A(1405) - 


-¥ Aj 


0.0011±0.0001 


1.65±0.031 


0.0223±0.002 


A(1405) - 


^£° 7 


0.170±0.012 


0.191±0.009 


0.437±0.013 


£(1385) - 


-»■ A7T 


1.421±0.0278 


0.0321±0.002 


0.0312±0.002 


£(1385) - 


-> £ + 7T- 


0.161±0.01 


0.00254±0.001 


0.00138±0.0006 


£(1385) - 


-> A7 


0.0184±0.002 


2.335±0.039 


0.0704±0.005 


£(1385) - 


-^£°7 


0.191±0.011 


0.058±0.0001 


0.225±0.015 


AK*+ -> 




0.213±0.010 


0.010±0.006 


2.931±0.051 


AK*+ -»■ 


iT+7 


0.0022±0.0001 


0.158±0.003 


2.351±0.046 



BACKGROUND CONTRIBUTIONS 

Tabic fll] lists all decay channels taken into considera- 
tion and the value of the acceptance for the confidence 
level cuts P%{x 2 ) < 1% followed by P*(x 2 ) > 10% 
for the 7- hypothesis and -P 7 (x 2 ) < 1% followed by 
P%(x 2 ) > 10% for the 7r°-hypothesis. To use these ac- 
ceptance terms to correct the signal yields, an estimate 
of the number n\ for the A(1405) in the event sample is 
required. The corrections for the 7 channel are given by 



N 7 (A* -> A 7 ) 
N 7 (A* -» E° 7 ) 
N y (A* -> SV) 
N 7 (A* ->■ E+tt") 



A$(A"f)BR(A* -y A 7 )n A 



AA 7r (S0 7r o )+A A 7r(s+7r - ) 
^(S° 7 )Bi?(A* -> £° 7 )n A 



ZT, (10) 



AA 7r (S0 7 rO)+AA 7r (S+ 7 r-) 
A A (£°7r> A 



zf(H) 



rr, (12) 



A^E+tT 



TT, (13) 



where BR is the branching ratio for the decay shown, 
and likewise for the ir° channel. 

Isospin symmetry is assumed so that BR(YPtt°) = 
BR(E+ir-) = BR{T,-tt+) w 1/3 for the A(1405) de- 
cay channels. The subscript "771" denotes the accep- 
tance for events that do not satisfy the confidence level 
cuts for either hypotheses of the kinematic fit {i.e. it 
is likely to come from some background reaction). The 
values for Bi?(A(1405) -> A 7 ) = 5.4 ± 0.2 x 1CT 4 % and 
BR(MU05) -> E° 7 ) = 2.0 ±0.1 x 10~ 4 % are taken from 
Ref. 0. 

Contributions from the E-7T decay of the E°(1385) are 
also considered. The term in Eq. (O that takes the E* — > 
E + 7r~ counts into consideration uses R^ — 0.135±0.011 
with the acceptance for the individual channels sub- 
ject to the radiative (tt°) hypothesis and denoted as 
A^Ett) (A^(Ett)). The E* -> E° 7 decay branch is also 
considered, however no measured branching ratio cur- 
rently exists for this channel. The acceptances are very 
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FIG. 11: Missing mass of jp — > K + A(X) for data (points 
with error bars) and Monte Carlo simulations for the 7p — > 
A' + A(1405) reaction (histogram) normalized to the data. 



small and using the higher theoretical prediction of the 
Algebraic model [20[ yielded only negligible contributions 
to the background. 

In order to find 7ia , one can look at the events for which 
neither the 7 nor the 7T° hypothesis is satisfied. The 
value of ha is difficult to determine due to the non-Breit- 
Wigner shape of the A(1405) decay. A better approach is 
to use a Monte Carlo simulation to fill the background ac- 
cording to its internal decay kinematics and normalize it 
to the data such that the MC matches the data, thereby 
giving an estimate of n\. Figure [TT] shows the MC sim- 
ulation normalized to the data, giving the estimate used 
for n\. This can be used to correct all backgrounds ex- 
cept for the K* . The final estimate found is n A — 4085. 

The 7p — > K*°T, + reaction was investigated with the 
MC simulation and compared with data. This back- 
ground was determined to have a negligible effect on the 
final result, since there is no A in the final state. For the 
■yp — > K* + A reaction, few events survive all of the cuts. 
To include corrections for the few events that do survive, 
an estimate of the K * + background must be established. 
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FIG. 12: Missing mass off the A fit with a Gaussian for the 
mass off of the K + in the window of 1.35-1.5 GeV. 



TABLE III: Breakdown of statistics for each term in Eqs. (8)- 
(9) for the A(7) and A(-7r°) hypotheses. Each listed channel 
is subtracted from the raw counts found from the kinematic 
fit in each case. The uncertainties are statistical only. The 
E(1385) — > E + 7T~ contributions are of the order of 10 -5 and 
smaller and are not included in the table. 



Reaction 




A( 7 ) 


A(7r°) 


Raw counts 


635 


13950 


A(1405) - 


-►eV 


3.41±0.36 


168.94±11.65 


A(1405) - 


^E+^" 


4.44±1.09 


98.98±7.80 


A(1405) - 


A7 


3.04±0.59 


O.OliO.OO 


A(1405) - 


^E° 7 


0.13±0.04 


0.12±0.04 


AK*+ -> 


K+7V° 


O.OOiO.OO 


27.41±1.71 


AA'*+ -> 


K+j 


0.03±0.00 


O.OOiO.OO 


Corrected counts 


623.92±25.23 


13654.56±118.95 



The correction for this background has the form 



RESULTS 



A„o{K 



*+ 



N T (K* K+tt") -- 
K+-K a )n(K*+ -)■ K+n") 



(14) 



A n (K*+ -)• K+n°) 



where A n (K* + — > K +/ k°) is the acceptance for the 
K* + — > K + tt° channel under the 7r°-hypothesis while 
A(K*+ -)• K+tt°) is the acceptance of the K*+ K+ir° 
channel which is dependent on the extraction method to 
obtain the K*+ -)• K+n counts. n(K*+ -> K+ir°) is 
the estimated number of K* + — > K + tt° events in the 
data sample. Similarly, the radiative decay of the K* 
has the form 



N 7 (K* ->■ Kj) = -R{K*+ -> K + i) 



(15) 



Ay(K* 



A{K*+ -> K+tt q ) 



K + l) 



with N W (K* Ktt°) from Eq. $T5§ and BR{K*+ -> 
K + j) ~ 9.9 x 10 -4 . An estimate of the number of K* 
events was obtained by matching the MC simulations to 
the data. The K* + — > K + tt° mass distribution has been 
fit as shown in Fig. [T^l In addition a fit is more easily 
obtained over the range that allows the higher part of the 
excited state mass spectrum to pass through. Using the 
resulting Gaussian fit of the K* peak while studying the 
mass off of the K + for mass windows ranging from 1.34- 
1.5 GeV to 1.34-1.8 GeV, we are able to extrapolate down 
to the nominal S*° mass cut (1.34-1.43 GeV), see Fig. g] 
Both methods gave similar results for n(K* + — ) K + tt°) 
used for the background correction in the final ratio. The 
extrapolated number of K* — > K + ir° events present were 
1207, 14.8 of which passed the 7r°-hypothesis. 

Each background contribution is enumerated in Table 
IIII The table gives the final breakdown of statistics for 
each term used in Eqs. (8)-(9) to achieve the corrected 
counts. 



The resulting corrected counts An 7 = 623.92 ± 25.23 
and An n = 13654.56 ± 118.95 are used in Eq. to 
obtain a ratio of 1.42 ±0.12(stat)%. The value of each 
cut was varied to study the effect on the final acceptance 
corrected ratio. For each variation the new acceptance 
terms along with the background contributions in Eq. 
(O were also recalculated. Each major systematic uncer- 
tainty contribution is listed in Table HVl as described in 
detail below. 

Particle identification was done by calculating the ve- 
locity for a particle of a given mass, using the measured 
momentum, and compared with that expected from the 
measured path length and time-of-flight. Kaons, protons 
and pions are selected based on the difference of the cal- 
culated and measured velocity, called A/3. Variations of 
the width of the A/3 cut to identify particles gave slightly 
different values of the final ratio, shown in line (1) of Ta- 
ble nvi 

The distance of closest approach (DOCA) cut for the 
proton and ir~ , used to reconstruct the A momentum, 
was varied and the stability of the final ratio was exam- 
ined. In the stable region, corresponding to cuts in the 
range from 3 cm to 14 cm, the effect on the final ratio 
is mostly toward higher values, for a larger DOCA cut 
value, as listed in line (2). 

Similarly, the value of the cut on the transverse mo- 
mentum, P xy , was varied. The ratio stabilizes starting 
at the cut point shown in Fig. [7J A series of cuts were 
used starting at 0.0009 GeV 2 and ending at 0.0025 GeV 2 
to study the effect on the final ratio, given by line (3). 

The Monte Carlo simulations for various background 
reactions were done assuming a i-dependent slope of 2.0 
GeV~ 2 for the differential cross sections, based on Regge 
theory, as described earlier. The value of the i-slope is 
not known precisely, and was varied by ±25%. The effect 
on the final ratio is shown in line (4). 
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The number of counts for the A(1405) and K* back- 
grounds were determined from fits to the data, using 
comparisons of shapes from Monte Carlo simulations 
with the data shown in Figs. llll and [T21 The uncertainty 
in the number of counts for these backgrounds also affects 
the final ratio, as shown in lines (5) and (6). 

To look at the systematic dependence on the choice of 
the confidence level cuts, the range defined by the uncer- 
tainty for P%(x 2 ) was checked. As previously described, 
the Monte Carlo studies lead to the set of optimal P" 
cuts for a given Pz. These optimal cuts allow one to re- 
cover, in our standard analysis framework, the ratio of 
A7 to Air decay that was used as input into the Monte 
Carlo simulations. To maximize counting statistics while 
minimize this uncertainty, the optimal cuts chosen for 
the analysis were < 1% and P% > 10%. 

The variation in the branching ratio is studied by se- 
lecting the confidence level cuts that lie slightly outside 
the optimization region found in simulations. In this way 
the largest range for P* and P° can be tested while still 
respecting the cuts derived from the optimization map. 
Using the full range of ratios in Table [V] the largest and 
smallest values show the variation for different choices of 
confidence level cuts, and is listed in line (7). 

The branching ratio of the radiative decay of the 
A(1405) affects the final result. Because this value is 
not measured directly, but is taken from the calculated 
value in Ref. [29j], there is some uncertainty in it. Using 
the range of values for this branching ratio given in the 
literature, and recalculating its effect on our result, leads 
to the uncertainty quoted in line (8). 

Table IIVI shows a summary of the systematic studies 
and the higher and lower value of the extracted ratio 
based on the variations mentioned for each type of un- 
certainty. The deviation of the ratio is defined by the 
difference from the quoted ratio of R = 1.42%. These 
deviations, shown in columns 3 and 5 of Table IIV| are 
added in quadrature to give the total systematic uncer- 
tainty. 

The range of the systematic uncertainty for R in Ta- 
bic IIVI is smaller than the statistical uncertainty, in part 
because each combination of cuts has a large overlap of 
events (i.e. the same subset of events is present for all 
choices of cuts). Since the kinematic fit requires a con- 
straint on the A mass, the kinematic cut on the invariant 
mass of the pn" has no effect. 

The final calculated ratio, given in percent, is 



R 



A 7 _ r[S°(1385)->A7] 
Att ~ r[£°(1385)-s-A7r"] 

= 1.42±0.12(stat)i°;^(sys)%. 



(16) 

Previously published work (llj on this branching ratio 
yielded a ratio of 1.53d=Q.39±g^f %. The value given here 
is consistent within uncertainties of the previous value, 
but has smaller uncertainties. The smaller uncertainty is 
important, as the previous uncertainty was on the same 



order as the theoretical meson cloud corrections to the 
EM decay of the A. If similar meson cloud corrections 
are to be proven true for EM decay of the S*° baryon, 
then the smaller experimental uncertainty is a significant 
improvement. 

The width for the branching ratio achieved comes 
from the use of the full width of the which is 

r(E*°)FuU = 36.0 ± 5.0 MeV with the branching ratio 
that the radiative signal is being normalized to, which is 
the P(£*° -> Att ) = 87.0 ± 1.5% [4]. The partial width 
calculation is then 



r s «o. 



>A7 - P^P(E*° -> A7r )r(E* ) Full 
445 ± 80 keV, 



(17) 



where a systematic uncertanty in R^ of ±0.11% is use in 
combination with the statistical uncertainty. Note that 
a large part of the uncertainty in Eq. (II T[) comes from 
the uncertainty of the full width r(E* )Fuii- 

These results verify that the partial width is indeed 
significantly larger than leading theoretical predictions, 
indicating that meson cloud effects are an important con- 
sideration for future calculations. 

The radiative decay E*° — > A7 is made up of Ml and 
E2 electromagnetic transitions. Assuming that the E2 
amplitude is very small, one can calculate the transition 
magnetic moment from the measured radiative width 
301, 



'2M p 2 r s , MAT 



ap^ 



fJ-N 



2.75 ±0.25/^, (18) 



where •p 1 is the photon momentum, M p is the mass of the 
proton, a = e 2 /47r ~ 1/137, and /xjv is the nuclear mag- 
neton. The value for the transition magnetic moment 
is larger than most model predictions even within the 
experimental uncertainty. For exam ple the naive quark 
model predicts ^s*°^A 7 = 2.28^jv UM- We hope that 
this measurement, along with others to come, will moti- 
vate theorists to understand the effect of the meson cloud 
on the magnetic moment and hence extend our knowl- 
edge of the quark wavef unctions in the decuplet baryons. 
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TABLE IV: Ranges of systematic variation of the final ratio, given in percent, along with the deviation from the central value. 



Source 


Low Value Low Deviation Hi 


gh Value Hij 


di Deviation 


(l)Particle identification 


1.380 


-0.040 


1.490 


+0.070 


(2)ptt~ DOCA cut point 


1.350 


-0.007 


1.480 


+0.060 


(3)Transverse momentum P xy 


1.415 


-0.005 


1.433 


+0.013 


(4)Monte Carlo t-dependence 


1.380 


-0.040 


1.440 


+0.020 


(5)A(1405) counts 


1.420 


-0.000 


1.470 


+0.050 


(6) K* counts 


1.420 


-0.000 


1.431 


+0.011 


(7)P(x 2 ) cut points 


1.388 


-0.032 


1.448 


+0.028 


(8)A(1405) — > A-y correction 


1.390 


-0.030 


1.420 


+0.000 


Total Uncertainty 




-0.072 




+0.112 



TABLE V: Dependence of corrected branching ratio for vari- 
ation of the confidence level cuts shown. 

Py(%) P°(%) R(%) 



15 


7.5 


1.388+ 0.12 


15 


5 


1.390+ 0.12 


10 


5 


1.422+ 0.12 


10 


1 


1.420+ 0.12 


10 


0.5 


1.421+ 0.12 


5 


0.1 


1.448+ 0.12 


5 


0.05 


1.436+ 0.12 
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Appendix: Ratio Derivation 

To calculate the ratio in Eq. (J7J, the leakage of the 
7T° region into the 7 region (and vice-versa) is the dom- 
inant correction. Taking just these two channels into 
consideration, the number of true counts can be repre- 
sented as N(Aj) for the E*° -+ A7 channel and N(Aw) 
for the E*° — »■ A-tt channel. The acceptance under the 
E*° — > A-f hypothesis can be written as A 7 (X), with the 
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subscript showing the kinematic fit hypothesis type and, 
in parentheses, the channel used in the Monte Carlo for 
the acceptance. For example, the calculated acceptance 
for the £*° — > A 7 channel under the £*° — > A 7 hypothe- 
sis is A 1 (Aj), whereas under the £*° —5- Att hypothesis 
it is A 7r (A'j). It is now possible to express the measured 
values n 7 and 

n 7 = A 7 (A 7 )7V(A 7 ) + A 7 (Air)N(Air) (19) 



n n = M^)N(An) + A ff (A 7 )JV(A 7 ). 



(20) 



The desired branching ratio of the radiative channel 
to the 7T° channel using the true counts is then R = 
N(Aj)/N(An). This can be obtained by dividing Eq. 
(fT9")) by Eq. (j2"0")) expressed in terms of R as 



n 7 _ iL4 7 (A 7 ) + A 7 (A7r) 
~~ A n (Air) + i?Ar(A 7 ) : 



(21) 



then solving for R. Expressed in terms of measured val- 
ues and acceptances, the branching ratio is 



R 



n 7 A T (A7r) — n,rA 7 (A7r) 
ra, r A 7 (A 7 ) - n 7 A 7r (A 7 ) 



(22) 



Equation (|22[) uses the assumption that contributions 
from the £(1385) will only show up as A 7 or Air , ne- 
glecting the £(1385) — > £7r channel. An estimate of the 
total number of £(1385)'s produced using the A7r° chan- 
nel is 



7V(£*°) 



A^° 



i?(£* ->• Att ),4(£* -> Att )' 



(23) 



where i?(£* — > An ) is the branching ratio of the £(1385) 
decay to An and A(£* — > Air ) is the acceptance for 
that channel. An estimate of the number of £(1385) — > 
£ + 7r~ — > pn°ir~ counts that would contribute to the 7r° 
peak is then: 



i?(£* 



N(T,* £^ 
£+tt-)A(£* - 



■R(S*-»S+x~)>l(£*- 
iJ(E*^A7r»)A(S*- 



►Att ) 



r") 

£+7r-)A r (£*°) 

iV(£* ->• Att°) , (24) 



where i?(£* — >• £ + 7r ) is the branching ratio of the 
£(1385) to decay into £+7r" and A(£* -)■ £+7r _ ) is the 
corresponding acceptance after all cuts. It is possible to 
simplify the expression by using, 



Rak - 



i?(£* 



£±7^1 



0.135 ±0.011, 



i?(£* -> Att ) 

using the PDG average value Q ■ The two charged com- 
binations of the £-7r decay have equal probability. The 
Clebsch-Gordon coefficient for the £* — > £ 7r° decay is 
zero, assuming isospin symmetry. The observed counts, 



expressed in terms of true counts and corresponding ac- 
ceptances for each hypothesis, becomes 



H A-n 



n 7 = A T (A 7 )AT(A 7 ) + (A 7 (Att) + -^-A 7 (Ztt))N(Att) 

(25) 



and 



(At(Att) 



H Ak 



A^tt))N{Att) + (^ r (A 7 ))JV(A 7 ). 

(26) 

Solving for R will result in a branching ratio that includes 
all needed information from the £(1385). Although the 
corrections to R from other contamination should be 
small, it is necessary to include them in the calculation. 
There is some probability that contamination for these 
other channels can leak through, and acceptance stud- 
ies were done for all channels under both the A 7 and 
An hypotheses. Results from the acceptance for each 
hypothesis are shown in Table [II] The branching ratio 
must include corrections for the K* + — > K + X and the 
A(1405) — > £-7r contamination, as well as a contribution 
to the numerator of R from A(1405) —> A 7 decay. The 
leakage of the £ 7 channel is assumed to be small relative 
to the A 7 signal. However, this channel is still considered 
in the acceptance studies, see Table UT1 

The branching ratio, taking these backgrounds into 
consideration, is Eq. §7$. The n 7 (n^) terms come di- 
rectly from the yield of the kinematic fits and represent 
the measured number of photon (pion) candidates. A 
similar notation is used so that the pion (photon) chan- 
nel identifications are denoted ^(£ + 7r~) (A^(Y; + tt~)), 
where A^(Att) is the relative leakage of the Att channel 
into the A 7 extraction, A^ (A 7 ) is the relative leakage of 
the A 7 channel into the Air extraction, A E is the accep- 
tance strictly for the £(1385), and A A is the acceptance 
for the A(1405). 

TableHTlshows the acceptance terms for the other back- 
ground channels that are considered in the ratio calcula- 
tion. The table lists three columns sorted by hypothesis: 
A-,, A n , and the counts that made all other cuts but did 
not satisfy either the 7 or ir° hypothesis, A 17! . The latter 
is used to obtain an estimate of counts for the specific 
backgrounds listed. 

To obtain the values of N(A*) from Eqs. (10)-(13) 
that must be subtracted from n 7 or n n , the number of 
observed counts of A* — > £°7r° are used with the accep- 
tances from the third column and the acceptance of the 
background channel of interest. As an example consider 
A^(A* — > £ + 7r~) under the ir hypothesis. The number of 
observed counts n(£°7r°) above the 7r° peak is given by 



N(A* 



n A 



R{A* £°7r )A A (£°7r )' 



(27) 



The notation n A here is shorthand for n(£°7r°), while 
R(A* -4 £°7T ) is the probability that the A(1405) will 
decay to £°7r° and A A (£°7r°) is the probability that 
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this decay channel will be observed after all the applied 
cuts. Isospin symmetry is assumed so that i?(E°7r°) = 
i?(S + vr-) = fl(£-7T+) w 1/3 for the A(1405) decay chan- 
nels. An estimate of the number of counts in the 7r° peak 
coming from the reaction A* — > £+7r~, using Eq. (1271) . 
is 



i?(A* -> E+7r-)A^(S+7r-)iV(A*) 
A$(X+n-)n A /A^(X°n°). 



(28) 



A small adjustment is made to ensure that the A* — > 
E + 7r~ contributions are also included by adding in 
the relative acceptance ^4 A (£ + 7r - ) to the denominator. 
These acceptance terms are found by independently using 



Monte Carlo for the -/p -> if +A(1405) ->• K+YPtt and 
7P ->• AT+A(1405) ->• if + S + 7T _ reactions. The counts 
that survive all cuts but did not satisfy cither the 7 or 
7T° hypothesis contribute to A 17V . The leakage for the 
7P — > _ftT + A(1405) — > AT+E+7T - channel is very small but 
is included for completeness. The final result is 



(E°tH>)+4* (E+7T-)- 



(29) 



From this example it becomes transparent how to express 
all other associated A(1405) corrections using only the 
observed n\ counts and the corresponding acceptance 
for this channel. 



